Non-additivity of decoherence rates in superconducting qubits 
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We show that the relaxation and decoherence rates T^^ and Tj"^ of a qubit coupled to several 
noise sources are in general not additive, i.e., that the total rates are not the sums of the rates due to 
each individual noise source. To demonstrate this, we calculate the relaxation and pure dephasing 
rates T^^ and T^^ of a superconducting (SC) flux qubit in the Born-Markov approximation in the 
presence of several circuit impedances Zi using network graph theory and determine their deviation 
from additivity (the mixing term). We find that there is no mixing term in and that the mixing 
terms in T^^ and Tj"^ can be positive or negative, leading to reduced or enhanced relaxation and 
decoherence times Ti and T^- The mixing term due to the circuit inductance L at the qubit transition 
frequency uiqi is generally of second order in ujoiL/Zi, but of third order if all impedances Zi axe 
pure resistances. We calculate Ti_2 for an example of a SC flux qubit coupled to two impedances. 



Introduction. The loss of quantum coherence and the 
transition from quantum to classical behavior has been a 
long-standing fundamental problem 1, 2J . More recently, 
the phenomenon of decoherence has attracted much in- 
terest in a new context, because quantum coherence is 
an essential prerequisite for quantum computation. For 
some systems that have been proposed as physical real- 
izations of quantum hardware (see, e.g., Ref. 0), there 
have been extensive studies, both in theory and experi- 
ment, of the mechanisms that are causing decoherence. 
Generally, an open quantum system loses coherence by 
interacting with a large number of external degrees of 
freedom (heat bath, environment). It is the physical na- 
ture of the environment and the system-environment cou- 
pling that distinguishes the various mechanisms of deco- 
herence. It is quite natural that for a given open quantum 
system there will be several distinct decoherence mech- 
anisms. Previous studies have typically tried to identify 
the strongest source of decoherence, i.e., the one that 
leads to the shortest relaxation and decoherence times, 
Ti and T2, and to analyze the corresponding mechanism 
in order to predict decoherence times. In the presence of 
several decoherence sources for the same system, the de- 
coherence rates T^^ and Tj"^ have usually been quoted 
separately for each source. Often, it is assumed that 
the total decoherence or relaxation rate is the sum of 
the rates corresponding to the various sources (see, e.g., 
Ref. for the case of superconducting qubits). In the 
theory of electron scattering in metals, this assumption 
is also known as Matthiessen's rule @. In this paper, we 
show that the total decoherence and relaxation rates of 
a quantum system in the presence of several decoherence 
sources are not necessarily the sums of the rates due to 
each of the mechanisms separately, and that the correc- 
tions to additivity (mixing terms) can have both signs. 

We investigate the decoherence due to several sources 
in superconducting (SC) flux qubits H S i, H E IH 
(see Ref. 4 for a review of SC qubits); the general idea 
of the present analysis may however be applied to other 
systems as well. SC flux qubits are small SC circuits that 



contain Josephson junctions. The differences ipi of the SC 
phases across the junctions Ji, where i = 1, . . . , n, are the 
relevant quantum degrees of freedom of the system; we 
denote the quantum operator of these phase differences 
collectively with the vector <p — (tpi, ip2, ■ ■ ■ , ^Pn)- The 
circuit is constructed such that it gives rise to a poten- 
tial U{(p) which forms a double well and therefore can be 
used to encode one qubit. In our analysis, we will make 
use of a recently developed circuit theory describing the 
dissipative dynamics of arbitrary SC flux qubits |l2l | . Our 
analysis relies on the theory for open quantum systems 
introduced by Caldeira and Leggett where the dissi- 
pative elements (impedances Zi) are represented by a set 
of baths of harmonic oscillators (an alternative approach 
to a quantum theory of dissipative electric circuits is to 
represent impedances as infinite transmission lines 'l3f|). 
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FIG. 1: Circuit graph of the gradiometer qubit [T3|, under 
the influence of noise from two sources Zi and Z2. Branches 
of the graph denote Josephson junctions Ji, inductances Li 
and Ki, current sources Isi, and external impedances Zi, and 
are connected by the nodes (black dots) of the graph. Inset: A 
resistively-shunted Josephson junction (RSJ) Ji, represented 
by a thick line in the circuit graph, is modeled by an ideal 
junction (cross) with critical current Id, shunt resistance Ri, 
and junction capacitance d. 
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For concreteness, we demonstrate our theory on the 
example of the gradiometer qubit with n = 3 junctions 
that is currently under experimental investigation |14| . 
see Fig. ^ We emphasize, however, that our findings are 
completely general and apply to arbitrary SC flux qubits. 
The qubit is controlled by applying a magnetic flux <i>c to 
the small loop on the left by driving a current Ibi in a coil 
next to it, and simultaneously by applying a magnetic 
flux $ on one side of the gradiometer using Ib2- Real 
current sources are not ideal, i.e., they are characterized 
by a finite frequency-dependent impedance ZijuA , giving 
rise to decoherence of the qubit 0, 0, 0, 0| . Since 
the shunt resistances Ri of the junctions are typically 
much larger (> Mf2) than the impedances of the current 
sources (between « 50il and « lOkil), we concentrate 
in our example on the impedances Zi and Z2 of the two 
current sources. 

Using circuit graph theory , we obtain the classical 
equations of motion of a general SC circuit in the form 



C03 = — — M * cp, 

dip 



(1) 



where C is the n x n capacitance matrix and 
U{ip; Ibi, IB2) is the potential. The dissipation matrix 
M(t) is a real, symmetric, and causal n x n matrix, i.e., 
M{t)^ = M{t) for all t, and M{t) = for t < 0. The 
convolution is defined as {f *g){t) — f{t — T)g{T)dT. 
Since it is not explicitly used here, we will not further 
specify U . The dissipation matrix in the Fourier repre- 
sentation [ill, M.{uj) = e-''^*-"*M(i) dt, where e > 
has been introduced to ensure convergence (at the end, 
e 0), can be found from circuit theory |ljj as 



-1 - T 

m , 



(2) 



where m denotes a real n x nz matrix that can be 
obtained from the circuit inductances, and where the 
nz X nz matrix (uj) has the form 



Lziuj) = Lziuj) 



(3) 



Here, nz is the number of impedances in the circuit (in 
our example, nz — 2) and Jjzi^^) = 'Z(uj)/iuj, where 
Z^uj) the impedance matrix. The frequency-independent 
and real inductance matrix Lc can be obtained from the 
circuit inductances ^2 ■ Since we start from independent 
impedances, Z and are diagonal. Moreover, note that 



ImL^ = Lu 



RcZ(tj) 



(ReZ(c^))-^Le(w) 



(4) 

where Lc(w) = Lc -t- lmZ{uj)/uj, thus it follows from 
ReZ > that ImL^^ and ImM arc positive matrices. 

Multi- dimensional Caldeira-Leggett model. We now 
construct a Caldeira-Leggett Hamiltonian |lj , H = Hs + 
'Hb+'Hsb, that reproduces the classical dissipative equa- 
tion of motion, Eq. and that is composed of parts for 
the system (S), for m > 1 harmonic oscillator baths (B), 



and for the system-bath (SB) coupling, 
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where the capacitor charges Q are the canonically con- 
jugate momenta corresponding to the Josephson fluxes 
(<i>o/27r)(p, where = (xai, ■ . ■ , Xam), and Ca is a real 
n X m matrix. From the classical equations of motion 
of the system and bath coordinates and by taking the 
Fourier transform, we obtain Eq. (Q, with M(ci;) 



(2V$o)2^^c„[m„(u;2_^2 



M((jj)-^, where the 



mxm mass and frequency matrices rcia and uJa are diag- 
onal with entries m^j and Uaj ■ Using the regularization 
UJ —^ u — ie when taking Fourier transforms also guaran- 
tees that M(t) is causal and real. 

Defining the spectral density of the environment as the 
matrix function 



where (X) — (5(Xy), we find the relation 
'$0 



3iu;) = 



27r 



ImM(ti)) = Jj{uj)m.j{uj)m.j{u 



(9) 



where we have used the spectral decomposition of the 
real, positive, and symmetric matrix 19] ImM(a;), with 
the eigenvalues Jj{uj) > and the real and normal- 
ized eigenvectors nij^Lu). The integer m < n,nz 
denotes the maximal rank of lmM(ci;), i.e., m = 
max^ (rank [ImM(Li;)]). Using Eq. and choosing 

Caij — JajT^ii^^aj), wc find that Jj{uj) IS the spectral 
density of the j-th bath of harmonic oscillators in the 



environment, Jj{uj) = (t^ /'i)Y.a{'^'ij I '^a]^aj)5{u 



The master equation of the reduced system den- 
sity matrix ps = Tr^p in the Born-Markov ap- 
proximation, expressed in the eigenbasis ||to)| 
of Tis, yields the Bloch-Redfield equation |2Cll |. 

Pnrn{t) = -iuJnrnPnm{t) ~ Y^kl ^nrnklPkl{t) , whcrC 

Pnm = (71 1 PS I m), a;„,„ = LOn - i^m, and cj„i is the 
eigenenergy of Tig corresponding to the eigenstate \m) 

The Redfield tensor has the form Rnmki — ^imj^r 
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J*_. = r|+)„fe: where HsBitUm = 

[n\e^*'^^HsBe '*'^^|to). For the system-bath inter- 
action Hamiltonian, Eq. (jT)), we obtain 
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FIG. 2: The relaxation rate Ti without the mixing term 
(dashed blue line), and including the mixing term for Rim — 
+ 10 kn (solid red line) and Rim = — lOkO (dot-dashed light 
blue line), for Afis = 0.5 pH as a function of KeZi. Inset: Ti 
for R = ReZi = 75 SI for a range of mutual inductances M13. 



where (^^^ — {n\ip\k). For two levels n = 0, 1, and within 
the secular approximation, we can determine the relax- 
ation and decoherence rates T^^ and T^^ in the Bloch 

equation as [12 Tf ^ = 2Re(r^|i^o + r^o(]i) and T^^ = 

(2Ti)-i + where T^' = Re(r(+^o + T[Vn ~ ^^Vii) 
is the pure dephasing rate. Using Eq. l(TU|l . we find 



Tfi = 4vl3MVoiCoth(^^y (11) 

= I lim>oo-<^ii)^^(¥'oo-¥'ii)- (12) 
With the spectral decomposition, Eq. @, we obtain 



Tfi = 4f]|v„i.m,(a;oi)|V,(a;oi)coth(^^), 



TZ' = |^K(0).(voo-<p,,)P^ 



(13) 
(14) 



In the last equation, we have used that the limit m.j (0) = 
limi^^om.j{Lu) exists because |mj(a;)p = 1 and thus all 
components of m.j{uj) are bounded. 

Mixing Terms. In the case where Lc is diagonal, or if 
its ofF-diagonal elements can be neglected because they 
are much smaller than lizii-^) for all frequencies uj, we 
find, using Eq. that the contributions due to dif- 
ferent impedances Zi are independent, thus m = nz and 
M(w) = mLz{i^)~^in^ = Y^^rhjiiv^iuj/{Zj{ijj)+iuLiLjj), 
where nij = rhj is simply the j-th column of the matrix 
m and Ljj is the j-th diagonal entry of Lc. As a con- 
sequence, the total rates l/Ti and l/T^ are the sums of 



1 a('^A \ - i2t3 ^oicoth (/3cjoi/2) 
-m = 4 — Iv'oi • mjl Re , (15) 
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1 



(16) 



In general, the situation is more complicated because cur- 
rent fluctuations due to different impedances are mixed 
by the presence of the circuit. In the regime Lc <^ lizi'-^), 
we can use Eq. (jSJ to make the expansion 



J-iz^ — i^z^ ^T-'z^'^c^'z^ ^ L^^'^LcL^^'^LcL^"'" — • 
The series Eq. (I17II can be partially resummed, 



(17) 



l,z\uj) = diag 



Zj{uj) + iujLj 



+ LmixH- (18) 



The first term in Eq. (|18|l simply gives rise to the sum of 
the individual rates, as in Eqs. (|15|l and H16() . while the 
second term gives rise to mixed terms in the total rates. 
The rates can therefore be decomposed as {X — 1,2,0) 



rpU) rT,(mix) 

j ^x 



(19) 



For the mixing term in the relaxation rate, we find 



1 



(mix) 



'^(^) '^oimImLj^(woi)m^¥'oiCoth^^^yi^ . 

(20) 

We can show that there is no mixing term in the pure 



dephasing rate, i.e., 

^(mix) nrp(niix' 



(mix) 



0, and consequently. 



T^"''-' = 2Tl^"""^ The absence of a mixing term in 
can be understood as follows. Since the first term in 
Eq. (|17|l only contributes to the first term in Eq. p8|l . 
the low-frequency asymptotic of ImLmix(w)~"'^ involves 
only o;^ and higher powers of u (it can be assumed 
that Zi{uj ~ 0) is finite), thus Eq. p2|) yields zero in 
the limit uj ^ 0. While ImL^^ is a positive matrix, 
ImL~?^ does not need to be positive, therefore the mix- 
ing term 1/T™"^ can be both positive or negative. Fur- 
thermore, we can show that if Z(w) is real, only odd 
powers of wLcZ"^ occur, and in particular, that in this 
case ImLniix(w)^^ = 0{uj^), by using Eq. |@J to write 
J(w) ~ wZ(w)-i - w3Z(u;)-iLcZ(w)-iLcZ(w)-i, up to 
higher orders in tt'LcZ(ct;)~^. 

In the case of two external impedances, nz — 2, we 
can completely resum Eq. (fT7|l . with the result 



4 



where are the matrix elements of Lc and where the 
approximation in Eq. (|21() holds up to 0{Z^^). In lowest 
order in l/Z^, we find, with ipi2 = (Voi ' ™i)(¥'oi ' ^2), 



1 



(mix) 



2 

Im 



Zi {uJoi)Z2{uJoi) 



coth 



(22) 



If Ri = ZilcuQi) are real (pure resistances) then, as 
predicted above, the imaginary part of the second-order 
term in Eq. ()21|) vanishes, and we resort to third order. 



mix jy TD I Lll I L22 L12 



Ri 



L22 

R.2 



(23) 



R2 



R2 



neglecting terms in 0{Rj^). If L12 <C Ljj, we obtain 
ImL-i^ « {u;^Li2/RiR2KLn/Ri+L22/R2)<J., and 



1 



(mix) 



271 J R1R2 \Ri R2)^ \ 2 

(24) 

For the gradiometer qubit (Fig. we find L12 ~ 
M12M13M34/L1L3, Lll ~ L2, L22 ~ ^4, where Lfc de- 
notes the self-inductance of branch {X=L or if) and 
Mfc; is the mutual inductance between branches Xk and 
X;, and where we assume Mij <^ L^. The ratio between 
the mixing the single-impedance contribution scales as 



(mix) 



Wqi-Li 



IIT\ 



(i) 



i?2 



i?, Lll 



L2 



(25) 



L, 



where we have assumed i?i i?2 

and (^01 ■ mi ~ V'oi ' ■ 

We have calculated Ti at temperature T <^ hujQi/kB 
for the circuit Fig. for a critical current Ic — 0.3 /lA 
for all junctions, and for the inductances Li = 30 pH, 
L3 = 680 pH, 2.2 = ^4 = 12 nH, M12 ~ n/TTL^, 



M34 ~ ^/L^Li (strong inductive coupling), — 6pII, 
with llIqi — 27r-30 GHz, and with the impedances Zi = R, 
Z2 = R + iRim, where R and i?ii„ = ilOkfi are real 
{Run > corresponds to an inductive, i?im < to a ca- 
pacitive character of Zi). In Fig.Ql we plot Ti with and 
without mixing for a fixed value of M13 = 0.5 pH and 
a range of i? = KeZi. In the inset of Fig. |21 we plot 
Ti (with mixing) and ((T^^^V^ + (Tf (without 
mixing) for i? = 75 51 for a range of mutual inductances 
Afi3; for this plot, we numerically computed the double 
minima of the potential U and c^qi for each value of M13. 
The plots (Fig. 12) clearly show that summing the deco- 
herence rates without taking into account mixing term 
can both underestimate or overestimate the relaxation 
rate l/Ti, leading to either an over- or underestimate of 
the relaxation and decoherence times Ti and T2. 

Higher-order terms in the Born series. Two series ex- 
pansions have been made in our analysis, (i) the Born 
approximation to lowest order in the parameter as ~ 
fiRg/Zi^LOQi) « l/wmTi, where /i is a dimensionless 
ratio of inductances and Rq = h/e^ is the quan- 
tum of resistance, and (ii) the expansion Eq. (|17|l in 
the parameter Q!l ~ UQiL/Zi, where L is the induc- 
tance of the circuit, where we included higher orders. 
The question arises whether the terms in the next order 
in as in the Born approximation could be of compara- 
ble magnitude to those taken into account in 1/tI™^^\ 
In our example, we could neglect such terms, because 
ctB/ctL ~ 0.001/0.1 — 0.01 1, but in cases where 
Q^B <^L^ higher orders in the Born approximation may 
have to be taken into account. 
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